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In this paper we present a theoretical and numerical analysis of the free expansion of a 
Bose-Einstein condensate, where we assume that the single particle energy spectrum is deformed 
due to a possible quantum structure of space time. Also we consider the presence of inter particle 
interactions in order to study more realistic and specific scenarios. The modified free velocity 
expansion of the condensate leads in a natural way to a modification of the uncertainty principle, 
which allows us to investigate some possible features of the Planck scale regime in low-energy 
earth-based experiments. 
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I. INTRODUCTION 

Over the last few years, the use of many-body systems as theoretical tools in the search for possible Planck scale 
effects has become a very interesting line of research HHS]. In particular, systems related to Bose-Einstein condensates 
are promising pathways in the search for low-energy traces of Planck scale physics [IHls]. These studies suggest that 
remarkable properties associated with Bose-Einstein condensates could be used to obtain representative bounds on 
the deformation parameters associated with quantum gravity models IBfllOj or in specific cases where it is possible 
to explore the sensitivity of these systems to Planck scale effects [MninHii]. 

An interesting phenomena related to Bose-Einstein condensates is the interference pattern when two condensates 
overlap [SHni- The interference pattern is a manifestation of the wave (quantum) nature of these many-body systems 
and could be produced even when the two condensates are initially completely decoupled. After switching off the 
corresponding traps, the system expand whilst overlapping, eventually producing interference patterns pRl - fT^ . 

Furthermore, it is well know that when the trapping potential is turned off the free velocity expansion of the cloud 
corresponds to the velocity predicted by Heisenberg’s uncertainty principle [HHizj. This fact is one of several reasons 
why Bose-Einstein condensates are relevant systems in the analysis and estimation of possible Planck scale effects, 
since quantum gravity models suggest modifications to this principle. |20H22] . 

Along these lines, in m we presented a study where the free velocity expansion of a Bose-Einstein condensate 
leads in a natural way to modifications of Heisenberg’s uncertainty principle. If we assume as a fundamental fact 
that the energy per particle is modified due to the quantum structure of space time, then the predicted modified free 
velocity expansion suggests a linear deformation in Heisenberg’s uncertainty principle, 

AxAp ^ ^ ~ + 0{a'^), (I) 

where a* = Mp is the Planck mass, c is the speed of light, and m is the mass of the particle. 

Additionally, is a real parameter of order unity, depending upon the quantum gravity model under consideration. 
As far as we know, this linear modification had not previously been reported in the literature, see for instance 
Refs. HOHSl]- 

The non-relativistic form of the aforementioned modified dispersion relation can be express in ordinary units as 
follows [23l I25j : 

The parameters ^1, ^^2; and ^3, are model dependent [131121], and should take positive or negative values close to 
I (see Ref. [13] for more details). In fact, the form of the energy dispersion relation ([^, was recently constrained 
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by using high precision atom-recoil frequency measurements |23l I25j . In this scenario, bounds for the deformation 
parameters of order ~ —1.8 ± 2.1 and —3.8 x 10® < ^2 < 1-5 x 10® were obtained. 

Eq.§ was deduced for a dilute system and neglecting the interactions among the particles within the condensate, 
i.e., the ideal case. Also, the modified Heisenberg’s uncertainty principle is a consequence of the leading order 
deformation contribution in Eq. ([^, which is linear in the momenta. In order to analyze the behaviour of the 
condensate under free expansion in the interacting case, more realistic scenarios are required. Clearly these must be 
taken into account when considering the corrections due by the next-to leading order term in Eq. (H. 

With this aim, we will analyze the behaviour of the solutions to this modified condensate scenario under free 
expansion using numerical tools, where we taken into account the effects produced by the leading order deformation, 
and the next-to leading order deformation in Eq. ([^ together with the interactions among the particles within the 
system extending the results reported in ca¬ 


ll. FREE VELOCITY EXPANSION OF THE CONDENSATE 


The modified energy associated with the system is given by 

Ei^if) = J dr 




(3) 


where is the wave function of the condensate or the so-called order parameter, V{r) = mujQr‘^/2 is the external 
potential, which is assumed to be an isotropic harmonic oscillator for simplicity. The term Uq = Airlifa/m, depicts 
the interatomic potential, a being the s-wave scattering length, i.e. only two-body interactions are taken into 
account. Notice that we have included in the total energy of the cloud, the leading order modification in the deformed 
dispersion relation Eq. M, through the linear operator |\/|v|®| [6l[T3j, where a = fimc/2Mp. Also, we have inserted 
the next-to leading order deformation in Eq. through the usual operator V, corresponding to the deformation 
parameter fi = ^ 2 / 2 Mp. Notice that this term is also quadratic in the momenta as is the corresponding kinetic energy. 
The corrections caused by the deformation parameter /?, could be re-absorhed in the usual kinetic energy term by 
defining the effective mass = Mpm/{Mp + f 2 m). We could also perform a similar analysis by assuming from 
the beginning that the deformation parameter fi is only a shift in the corresponding particle mass. However, as was 
pointed out in Ref. [B], both approaches lead to the same predictions for the ground-state energy and its properties, 
at least to first order in ^ 2 - Thus, without loss of generality, we analyze in this work the modifications caused by fi 
as independent contributions to the total energy of the system. 

Furthermore, we have assumed that ^3 = 0. Some insights about this latter parameter in the case when is non-zero 
will be given at the end of this letter. If we set a = /3 = 0 in the total energy Eq. ® we recover the usual expression 
[15]. 

The total energy of the cloud can be expressed as follows: 


E = Ep + Ep, (4) 

where Ep is the kinetic energy associated with particle currents 

= (5) 

The function Ep can be rewritten using the following components 

Ep = Eq + Ep + Ej -\- Ea + Ep, (6) 

where 



( 7 ) 
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so, {Eq) is related to the contributions of the ground state energy, (Ep) the contributions of the trapping potential, 
and (Ei) the contributions due to he particle interactions within the condensate. The contributions {Ea) and (iJ/j) 
contains the contributions of the deformation parameters a and /3, respectively. 

Firstly, Ep can be written as 


Ep = - 


3 


4 


N + -mojo^R^N + 


Uo 2 h ^ 3^2 

- - a—^N + B — -N, 

2{2T:f^^R3 V^R 2ii2 


( 8 ) 


where we have used the ansatz 

iVl /2 

= t: 3 / 4 r 3/2 exp(-r /2R ) exp[i(j){r)], (9) 

where N is the corresponding number of particles and ^(r) is a phase related to particle flows in the system. 

The choice of the ansatz ([^ , for the case of a weakly interacting Bose-Einstein condensate trapped in an isotropic 
three-dimensional harmonic-oscillator potential, seems to be a good conjecture for several reasons: First, Eq. ([^ 
clearly reflects the symmetry of the trap and in the non-interacting case is the exact solution of the corresponding 
equation of motion; Secondly, as was proven in the experiment described in the system operates deeper in the 
linear regime for sufficiently large expansion times, i.e., the system evolves almost as in the non-interacting case in 
this situation. This fact further supports the use of the ansatz In other words, in the experiment m was shown 
that the free velocity expansion at large times confirms that the evolution of the condensate can be independent of 
interactions during extended free fall experiments. Accordingly, the free velocity expansion can be computed in this 
scenario without loss of generality, by using the aforementioned ansatz at least to first order approximation in the 
deformation parameters a and /3. Thirdly, as we will show later on this paper, all these facts indicate that large 
expansion times are also relevant in the search for some Planck scale effects. 

Let us add that possible contributions due to the deformation terms can appear in the order parameter Eq. (1^ 
through the corresponding phase (j>{r), which is related to the local velocity of the condensate as v = [h/rnfV(j){r) 11^ . 
However, within our approach, all the measurable quantities of interest are computed by taking the norm of Eq. ® , 
which is related to the density and its derivatives (see Eqs.Q). In other words, it is necessary to calculate the full 
solution of the equation of motion, e.g., the corresponding Gross-Pitaevskii equation, together with the contributions 
due to the deformation terms. This general version of Eq. (|^ can be helpfully to analyze the contributions of the 
deformation terms upon the phase, the density and its derivatives. In order to test the validity of our model, the 
eventual predictions from the full solution can be useful to compare with the approach presented in this work. This 
is a non-trivial topic that deserves deeper analysis and it will be presented elsewhere. 

Let us start with our model by considering that the external potential V(r) is turned off at t = 0, in such a case 
there is a force that changes R and produces an expansion of the cloud [15]. It is straightforward to obtain the kinetic 
energy Ep by using the ansatz Eq. ([^, with the result Ep = “iRfNm/A. Moreover, assuming that the energy is 
conserved at any time, we obtain the following energy conservation condition associated with our system 


3TOi?2 3^2 UqN 2h 3^2 _ 3^2 U^N 2h 

4 4TOi?2 2(27r)^/^i?3 °^^/7^R~^ 2R'^ ^rnRl^ 2{2'Kf/‘^Rl y/f^Ro^ ’ 

where the dot stands for derivative with respect to time and i?o is the radius of the condensate at time t = 0, 
which is approximately equal to the oscillator length Oho = in the non-interacting case. Otherwise, when 

interactions are present, we will assume that the initial radius corresponds to the result for an isotropic trap [15j 


i?o — 



( 11 ) 


Additionally, R is function of time and corresponds to the radius at time t. If we set a = jB = Q then we recover 
the usual solution in the non interacting case |15j which is given by 




( 12 ) 


Notice that in the usual case, a = /3 = 0, vq = h/mRo, is defined as the velocity expansion of the condensate, 
corresponding to the velocity predicted by Heisenberg’s uncertainty principle for a particle confined within a distance 
Ro [I5|. Thus, in the usual case a = fB = 0, the width of the cloud at time t can be written in its usual form 

R^{t)=Rl + ivot)^. 


(13) 
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It is noteworthy to mention that when interactions are neglected we are able to obtain an analytical solution for 
Eq. (10 1 when R >> Rq together with a ^ 1 and /3 <C 1. In such s scenario we obtain 


^a,/3 W — ^0 + 




rrRRQ 


1 + 2m/3 — a 


8 h 

Si/tt mRo 




(14) 


If we set (3 = Q, the result obtained in Ref. m is recovered. Thus, we may recognize the free velocity expansion in 
function of the deformation parameters a and /3, which is given by 


= 




m^i?Q 


1 + 2mB — a 


ttiRq 

Since the corrections caused by a and j3 are quite small the following expansion is justified: 




h 4 

-(1 + 2m/3) — a- 


mRo'' ' ''' Sy/n 

Then, the velocity expansion corresponds to the following deformed Heisenberg’s uncertainty principle 

AxAp > ^(1 + /?*) — a*x + ..., 


(15) 


(16) 


(17) 


where we have defined 13* = ^ 2 'm/Mp and a* = ^i2TO^c/3v^Mp together with i?o = x. It is not surprising that the 


functional form of Eq. ( |17[ ) implies the following minimum measurable momentum and maximum measurable length 

Ap > (Apmin) « -a*, (18) 


, , h / I + p* \ 

Ax < (AXmax) « 2 ( Q,* j • 


(19) 


Notice that the inequality (18 1 is relevant only when < 0, which implies negative values of a*. These conditions also 
set the value range of deformation parameters and ^2 without breaking the inequality 0 - From a phenomenological 
point of view, these conditions can be used in other systems, in order to explore some issues related with the quantum 
structure of space time. 


III. NUMERICAL ANALYSIS 


In order to explore the velocity expansion and the possible corrections caused by the deformation parameters 
and ^ 2 , we need to solve Eq. (10 1 numerically at any time and taking into account the interactions among the 


particles within the system. We will take fiducial laboratory conditions over the parameters related to the model as, 
N ^ 10"^ — 10® particles, wq ~ 10 — 10® Hz, a ~ 10“® m, and m ^ 10“^® kg [55]. Additionally, Mp ~ 2.18 x 10“®kg, 
h ~ 6.623x 10“®'*Js, c ~ 3x 10®m/s. The deformation parameters considered here will be of the order of ~ —1.8±2.1 
and -3.8 x 10® < 6 < 1-5 x 10®. 

To illustrate how all these ingredients work together properly, we study now the accuracy of our numerical solution 
in some particular cases. These cases are: [I] For a,P and Uq non-zero. [H] For a = P = 0 and Uq = A-KhPa/m and 
[HI] For a, /3, Uq = 0. Also, we impose initial conditions for these cases. The condition at t = 0 for the cases [I] and 
[H] is described by the Eq.(ll |. The condition for the Case [HI] is set by R{t = 0) = {h/mwoy^^. In Figurej^we show 
the numerical solutions for R(t). Cases [I] and [H], (red and blue lines, respectively) show an identical evolution. The 
fact that the numerical solution for the Case [I] looks cut, shows that the code remains convergent a large times. Case 
(HI) represented by the green dashed-dotted line gives the exact solution R{t) = . 

The solution for this latter is illustrated in the the plot inside Figure 1. 

In Figure are illustrated the numerical solutions for the modified velocity R'(t) for the cases [I] and [H]. The 
velocity in Case [HI] is linear in time and its solution is exact, which for our interest in this figure with only show the 
modified cases. 

We notice interesting points related to the modified velocity and the computation of modified Heisenberg’s 
uncertainty principle in each scenario: 


• At early times, Eq. (171 shows a value around AxAp > 1.96 x 10“®®. This range correspond to the point where 
the numerical solutions of R'{t) for the cases [I] and [H] overlap as we see from the left plot inside Figure]^ 
This overlapping can be due to the equal dominance of both deformed parameters and ^2 ■ 
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FIG. 1: Case [I]: For a,/3 and Uo non-zero is represented by the red solid line. Case [II]: For a = /3 = 0 and Uo = ‘iTrh^afm 
is represented by the blue dashed line. Case [III]: For a,j3,UQ = 0 is represented by the green dashed-dotted line. The plot 
below shows the late time solution for Case [III]. 
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FIG. 2: Left-. For R'{t) at early times only the cases [I] and [II] are shown due the exact solution for Case [III]. Being so, Case 
[I] correspond to the red solid line and Case [II] correspond to the blue dashed line. The plot below represents the evolution at 
early times. Right-. For R'{t) at late times t > 4 sec, the plot below represents the evolution of R{t) in the same time range. 




























6 


• At the left of Figure.we observed a bounce in the solutions at early times. During this stage we have for Case 
[I]: AxAp > 1.28 x 10“^^ and for Case [II]: AxAp > 4.38 x 10“^®. 

• After this instability, the evolution of the modified velocity shows for Case [I]: AxAp > 4.59 x 10“^"^ and for 
Case [II]: AxAp > 7.22 x 10“^®. We observed a linear evolution of the modified velocity in where the modified 
Heisenberg’s uncertainty principle for Case [I] is larger than Case [II] due the eventual dominance of the deformed 
parameters and ^ 2 - This results was expected due the appearance of these deformations at small scales. 


At late times, let us consider, e.g. t « 4 sec, as in free fall experiments m- In this scenario we obtain for large 
times that AxAp > 10“®° (see at the right of Figure.]^. The evolution of the modified velocity for Case [I] is a 
constant. This behaviour is in agreement with the theoretical result that at large times the corrections caused 
by the deformation parameters could be representative and consequently, can be described by the modified free 
velocity Eq. (161 when interactions are neglected [13]. 


Generally, we notice that at very early times of expansion there is a period in which the velocity seems to be 
dominated by the deformation parameters. However, we estimate that this short period of expansion of order t = 
7.744 X 10“® sec may be hardly accessible from the experimental point of view, since some results offer order of 
milliseconds HH, i.e., three orders of magnitude bigger than the expansion time obtained here. Conversely, for large 
expansion times up to t = 4 sec, there is a region in which the presence of the deformation parameters modified the 
velocity expansion in a way that may be significant, even when interactions are present. 

Concerning the experiment performed in m, it was proven that for sufficiently large expansion times, the system 
operates deeper in the linear regime, i.e., in the non-interacting case. In consequence, the observed spatial interference 
pattern indicates that the fringe spacing scales linearly with the time of expansion and is inversely proportional to the 
initial separation of two condensates. In this experiment it was shown that the free velocity expansion at large times, 
confirms that the evolution of the condensate can be independent of interactions during extended free fall experiments. 
Each of the above scenarios shows that the modified free velocity expansion leads to deformations of Heisenberg’s 
uncertainty principle which are around two orders of magnitude smaller than the typical case. This fact, could be 
tested, in principle, in the laboratory, if the possible corrections in the free velocity expansion of the condensate can 
be eventually be measured. However, let us remark that according to our results, large expansion times are required. 
This analysis opens a very important branch of research concerning the search of some quantum gravity traces in low 
energy earth based experiments. 


IV. CONCLUSIONS 

We have analyzed and described the free velocity expansion of a Bose-Einstein condensate at different times and also 
when interactions are present assuming a deformed dispersion relation as a fundamental fact. Additionally, we have 
obtained a deformation of Heisenberg’s uncertainty principle which appears naturally just by looking up the modified 
free velocity expansion. However, for a further insight into this deformation the third deformation term in Eq. ([^ 
must also be taken into account. According to m this cubic term could be interpreted as inversely proportional to 
the lifetime of the condensate if we assume that this contribution to the total energy is imaginary. These facts, lead 
us to think that some of the particles leave the system (the condensate). In other words, this last assertion suggests 
that some of the particles forming the condensate, may be transferred to the excited states and in consequence could 
lead to instabilities within the system at some given time. Moreover, this deformation would contribute also to the 
functional form of the deformed uncertainty principle. These are non-trivial topics which deserve deeper investigation 
and on which we will report elsewhere. 

Finally, according to our results there are two relevant scales of time associated with the free expansion, which 
offers a possibility to detect small signals or traces from the quantum structure of space time. However, we stress that 
an optimal scenario in searching these possible signals is when the system expands for large times. As we mentioned, 
free fall experiments could provide signs of Planck scale physics in this scenario. 
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